Preprint typeset in JHEP style - HYPER VERSION Cavendish-HEP-10/17 

CERN-PH-TH/20 10-236 
DAMTP-2010-82 
MCnet/10/19 



Searching for third-generation composite 
leptoquarks at the LHC 



Ben Gripaios 1 , Andreas Papaefstathiou 2 , Kazuki Sakurai 2,3 , Bryan Webber 2 

1 CERN PH- TH, Geneva 23, 1211 Switzerland 
2 Cavendish Laboratory, University of Cambridge, Cambridge, UK 
^Department of Applied Mathematics and Theoretical Physics, University of 
Cambridge, UK 

Email: |gripaios@cern . ch , andreasQhep . phy . cam .ac.uk, 



sakuraiOhep . phy . cam . ac . uk , webber@hep . phy . cam .ac.uk 



Abstract: Fermion masses may arise via mixing of elementary fermions with com- 
posite fermions of a strong sector in scenarios of strongly-coupled electroweak sym- 
metry breaking. The strong sector may contain leptoquark states with masses as 
light as several hundred GeV. In the present study we focus on the scalar modes 
of such leptoquarks since their bosonic couplings are determined completely and 
hence their production cross sections only depend on their masses. We study all the 
possible gauge-invariant non-derivative and single-derivative couplings of the scalar 
leptoquarks to the quarks and leptons, which turn out to be, predominantly, of the 
third generation. We examine their phenomenology and outline search strategies for 
their dominant decay modes at the LHC. 
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1. Introduction 

Any model of electroweak symmetry breaking must explain not only the masses of 
gauge bosons, but also those of the Standard Model (SM) fermions. In models where 
electroweak symmetry breaking is driven by strong dynamics, it is usual to assume 
that fermion masses arise in much the same way as in the SM. That is, fermion masses 
arise via a Yukawa-type coupling, but in which the SM Higgs is replaced by some 
scalar operator of the strong sector, O, carrying the same gauge quantum numbers 
as the Higgs. However, this assumption necessarily makes it rather difficult, if not 
impossible, to satisfy constraints coming from flavour-physics experiments, whilst 
simultaneously providing a natural explanation of the hierarchy between the weak 
and Planck scales. The tension arises in the following way: suppose that flavour 
physics arises at some scale Af. In order to suppress dangerous flavour-changing 
processes coming from four- fermion operators, one needs Af > 10 3 TeV. The mass 
of the top quark arises from an operator of the form 

r r- OQ tC 

where d > 1 is the dimension of the operator O. To generate the large top mass 
with Af > 10 3 TeV requires d < 1.2 — 1.3 [1]. Now, for the hierarchy problem to be 
solved, the theory should not contain any relevant operators that can be added to 
the Lagrangian. In particular, the dimension of any gauge-singlet operator arising in 
the operator product expansion of O^O should roughly exceed four. This is clearly 
impossible in theories, such as weakly-coupled theories or strongly-coupled theories 
of large-N type, where correlation functions of operator products approximately fac- 
torize. Moreover, there are strong indications that it is also impossible in any theory 
which possesses a large hierarchy of scales and is thus is approximately a conformal 
theory in between [2-4]. 

An alternative, proposed long-ago [5], is to declare that the observed fermion 
masses arise by mixing of elementary fermions with composite, fermionic resonances 
of the strong sector. The latter are, of course, sensitive to the electroweak-breaking 
dynamics of the strong sector and transmit it to the elementary fermions via the 
mixing. Such a mechanism gives an automatic suppression of flavour- changing pro- 
cesses since the light fermions (for which the constraints from flavour experiments 
are strongest) are those which are least mixed with the flavour- changing dynamics 
of the strong sector. Moreover, it offers the hope that the observed hierarchies of 
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masses and mixings of SM fermions may be related to the electroweak hierarchy, via 
strongly-coupling effects. 

If fermion masses do arise in this way, it follows that the strongly-coupled sector 
knows not only about the SU(2) L x U(l)y of electroweak symmetry of the SM (as it 
must, since it breaks it to the electromagnetic U(l) em ), but also about the SU(3) of 
the colour interactions [6]. More precisely, the strongly-coupled sector must, at the 
very least, contain colour-triplet fermionic resonances that can mix with elementary 
colour-triplets to make the observed quarks. 

It is reasonable to expect that the strongly-coupled sector will contain other 
coloured resonances. If it contains bosonic, coloured resonances, then, depending on 
the particular gauge charges (or other global charges) these may be able to couple 
to a lepton and a quark, thus playing the role of leptoquark states. 1 For example, 
in the original model of [5], the fermionic resonances arose as technibaryons of a 
technicolour SU (3) interaction and (though not mentioned in [5]) such a model would 
contain leptoquarks arising as techni-mesons. 

Such resonances will typically lie around the TeV scale, but may even be sub- 
stantially lighter if they arise in the theory as pseudo-Nambu Goldstone bosons. 
Typically, such light states would be ruled out by flavour experiments, but in these 
models, leptoquark mediated flavour-changing processes are suppressed in exactly 
the same way as other flavour-changing processes that occur in the theory. Indeed, 
estimates given in [6] suggest that such states may, as a result of this flavour sup- 
pression, be as light as the current Tevatron bounds of roughly 200 GeV. 

These composite leptoquarks therefore make an ideal target for LHC searches, 
and if discovered would give us a strong hint about the mechanism of electroweak 
symmetry breaking. 

Since the leptoquark couplings to light fermions are highly suppressed, the only 
relevant couplings for direct collider production and detection are those to third 
generation fermions. 2 As a result, the leptoquark states will decay exclusively to 
third-generation fermions, that is to tr or tv T or br or bv T . Naively, since the lep- 
toquark couplings scale roughly with the Yukawa couplings, and since the Tevatron 
bounds preclude a leptoquark mass below m t , 3 one might conclude that decays in- 
volving the top must dominate. However, we shall see later that the gauge quantum 
numbers sometimes preclude couplings to top quarks and, of course, unknown global 
symmetries may also preclude one or more couplings. Thus we consider all four 
possible couplings. 

Since leptoquarks couple dominantly to third generation quarks and leptons, 
pair-production through colour gauge interactions will overwhelmingly dominate sin- 

1 Di-quark states are also possible, but we do not address them here. 
2 For an alternative scenario with leptoquarks of this type, see [7]. 

3 Searches at DO for third generation scalar leptoquarks decaying exclusively to br or bv T yield 
bounds of 210 GeV [8] and 229 GeV [9] respectively. 
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gle production at the LHC The channels of interest therefore involve pair-wise com- 
binations of tr or tv T or br or bu T . A The 262r and 2b + E(t channels already have 
been the subject of searches at the Tevatron [8,9], and can be adapted easily for the 
LHC. The use of novel kinematic variables such as My 2 in this 2b + EJt channel may 
well improve the prospects for discovery and mass measurement. The two channels 
involving the top require more ingenuity, but merit investigation. 

In the present paper we perform the first detailed phenomenological study of 
the possible production of such states at the LHC. In Section ^| we briefly review 
their quantum numbers, couplings and decay modes, which we have implemented in 
the general-purpose event generator HERWIG++ [10, ll]. 5 This allows us to propose 
and investigate some strategies for reconstructing third-generation leptoquark masses 
from their decay products, including those that involve top quarks, in Section |[ Our 
conclusions are presented in Section £|. 

Although we focus here on direct searches at the LHC, there are also promising 
channels for indirect searches, namely in — > K]x^ and B s — > fifi at LHCb, in 
// — > cy and r — > fi'j, in /i — e conversion in nuclei, and in r — > rj/j, at future B 
factories [6]. 

2. Phenomenology 

2.1 Scalar leptoquark pair-production 

We focus on scalar leptoquarks in the present study since their bosonic couplings 
are determined completely by QCD and hence their production cross sections only 
depend on their masses. Moreover, the lightest (and most easily accessible) lep- 
toquarks in these scenarios arise as scalar, pseudo-Nambu Goldstone bosons. The 
type of leptoquarks we are considering are predominantly pair-produced via gluon- 
gluon fusion or quark-antiquark annihilation, due to the fact that they couple to the 
third generation quarks and leptons. Only charge-conjugate leptoquarks can be pro- 
duced in this way: associated production of different leptoquarks is forbidden since it 
would not conserve the Standard Model gauge quantum numbers. Single production 
in association with a lepton is allowed but at a 14 TeV LHC it becomes dominant 
at leptoquark masses of about 2.2 TeV, at which point the total cross section is 
a ~ 10~ 2 fb, already too low for discovery. 

2.1.1 Effective Lagrangian for interactions with gluons 

The effective Lagrangian describing the interaction of the scalar leptoquarks with 
gluons is [15] 

C? s = (Dtj& y (Iff $») - M^tfr. j (2.x) 

4 Note that third-generation lepton-quark couplings are also possible in i?-parity- violating super- 
symmetric theories. 

5 PYTHIA [12-14] contains an implementation of a single scalar leptoquark of arbitrary flavour. 
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where $ is a scalar leptoquark, i,j, k are colour indices, the field strength tensor of 
the gluon field is given by 

Q% = d,A a v - d v Al + g s f abc A, b A uc , (2.2) 

and the covariant derivative is 

= - igJ%Al • (2.3) 

The Feynman rules that result from this Lagrangian and the diagrams that contribute 
to pair-production of scalar leptoquarks are given in Appendix [A]. Expressions for 
the cross sections are given in Appendix |B[ 

2.2 Leptoquark decays 

2.2.1 Non-derivative fermion couplings 

The effective Lagrangian that describes the possible non-derivative couplings of the 
scalar leptoquarks to third-generation quarks and leptons is given by [16]: 

£nd = (goLol^h + goRt c R T R )So 

+ (h 1L u R £ L + h 1R q L iT 2 T R )S 1/2 + h 2L d R £ L S 1 / 2 + h.c. , (2.4) 

where the r a are the Pauli matrices, gt and £l are SU(2) L quark and lepton dou- 
blets respectively and t R , b R and r R are the corresponding singlet fields. We denote 
charge conjugate fields by f RL = (-Pr,l/) c , where the superscript c implies charge 
conjugation. In Table [I] we give the quantum numbers for the five types of non- 
derivatively coupled scalar leptoquarks: the 5'L r (2) i -singlet complex scalars 5*0, 5o, 
the 5'[/(2) L -triplet complex scalar Si and the S"?7(2) L -doublets 5*1/2 and Si/ 2 . 

The numbering scheme used in our implementation of scalar leptoquarks in the 
HERWIG++ event generator is also given in Table |]. The particles are numbered as 
99NDDDJ, where iV distinguishes the representation of the standard model gauge 
group, DDD is the lowest possible number chosen to relate the leptoquark to the 
Particle Data Group (PDG) codes of decaying fermions, and J = 2S + 1, where 5* is 
the particle spin. The sign of the PDG code is negative for colour anti-triplets and 
positive for colour triplets. Hence —9911561 is the 'first' type of leptoquark, 5*0, and 
can decay to particles with codes 15 (r) and 6 (t). 

Notice that the first three kinds of leptoquarks, the So, Sq and the Si triplet are 
colour anti-triplets and the particles (as opposed to the anti-particles) decay into an 
anti- lepton and an anti-quark. This is contrast to the S1/2 and S1/2 doublets, which 
are colour-triplets and decay into quarks and anti-leptons. 
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Table 1: Numbering scheme, charges and possible decay modes for the non-derivatively 
coupled scalar leptoquarks. Y represents the U(l)y charge and T 3 is the third component 
of the SU(2)l charge. Since Si is an SU{2)l triplet, it contains three complex scalars. The 
S1/2 and Si/2 are SU(2)l doublets. The naming convention is explained in the text. The 
minus sign in the ids of some of the leptoquarks indicates the fact that they are anti-triplets 
of SU(3) C . 

2.2.2 Derivative fermion couplings 

We also consider leptoquarks that couple derivatively to the quarks and leptons. 
The couplings of the leptoquarks to fermions involve three fields, and hence two 
independent positions for the derivative to act, modulo integration by parts. Here, 
we choose to put the derivative on either the quark or the lepton, such that the 
Lagrangian is given by: 
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(2.5) 



where the index a G {1,2,3} and p^ 1 , i 6 {I, q}, denotes the momentum of the lepton 
or quark. 

The charges of the primed scalar states appear in Table Q they correspond, 
of course, to those of vector leptoquarks. Notice that whereas the So is a colour 
anti-triplet, S' is a colour triplet and so on. 

Consider a leptoquark S' that couples derivatively to fermions in the following 
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Table 2: Numbering scheme, charges and possible decay modes for the derivatively- 
coupled scalar leptoquarks. The details are as in Table 



way: 
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(2.6) 



where the / is the sigma-model scale for the strong dynamics. Consider the decay 
of the S' Q to on-shell fermions via the coupling g' QL JjlP^I^l- We then have: 



QolM^'o = 9' L, q bLf S oh + Sol, A A^o 

= 9oL,q m bbRSo£L + g' 0L/ m T b L S' r R 



(2.7) 



Note that the chirality of one decay product is reversed in each term by the mass 
insertion, which breaks the gauge symmetry. An equivalent manipulation is given in 
Appendix for terms that contain conjugate fields. For simplicity of the analysis, we 
choose to set the quark and lepton primed couplings for each term equal, g' e = g' = g', 
where g' represents g' , g[ or h[, for the rest of the paper. As a result of the above 
manipulation, an effective Lagrangian for the on-shell decay of a scalar leptoquark 
S' Q may be written as: 
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Table 3: The Aj couplings of the non-derivatively scalar leptoquarks to the different 
quark-lepton combinations, as they appear in the Lagrangian. 

effectively converting all the derivative couplings to ones that look like those for the 
unprimed leptoquarks, with the lepton or fermion masses appearing in the coupling. 
See Appendix [C] for the full effective Lagrangian. Since the scale / is typically 
a few hundred GeV, couplings proportional to the top quark mass are expected 
to dominate when the corresponding decays are kinematically allowed. The on- 
shell fermion assumption is realistic since the widths of the fermions are small in 
comparison to their masses and hence off-shell effects are negligible. 

2.2.3 Decay widths 

The decay width of non-derivatively coupled scalar leptoquarks in the limit of mass- 
less quarks and leptons can be calculated by [16]: 

F = ^ (Xl(tq) + Xl(*q) + X&vq)) , (2.10) 

where the couplings \l,r(£q) for the types of leptoquarks we are considering are given 
in Table fj in terms of the couplings that appear in the Lagrangian. The couplings 
are taken to be real. The expression gives, for quark-lepton couplings g ~ 0.1 and 
leptoquark mass of ~ 400 GeV, a width of ~ 0.1 GeV. The decay width to massive 
q£ is further suppressed by a phase space factor compared to the massless quark and 
lepton width [9]: 

F ~ (1 - r q - r e )^l + (r q -n) 2 -2r q -2n , (2.11) 

where r q / are the squared ratios m^/M| Q respectively. 

Table ^ shows the couplings for the primed, derivatively-coupled, scalar lepto- 
quarks. The expression for the width given in Eq. ( |2.10| ) remains unchanged in the 
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Table 4: The Aj couplings of the derivatively-coupled (primed) scalar leptoquarks to the 
different quark-lepton combinations, as they appear in the Lagrangian. In our analysis, we 
have set the quark and lepton couplings equal for simplicity. 

case of the primed leptoquarks, with the couplings Aj taking the appropriate values. 
Tables |5] and ^| show example decay widths and branching ratios for scalar lepto- 
quarks of mass Mlq = 400 GeV and couplings g = 0.1. In the case of derivatively 
coupled leptoquarks we choose a suppression scale / = 800 GeV. 
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Table 5: Decay widths for non-derivatively coupled scalar leptoquarks of mass Mlq = 
400 GeV and couplings g = 0.1. 
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Table 6: Decay widths for derivatively-coupled (primed) scalar leptoquarks of mass 
Mlq = 400 GeV, couplings g' = 0.1 and suppression scale / = 800 GeV. 



3. Experimental search strategies 
3.1 Reconstruction strategies 

Table [7| provides an overview of our suggested reconstruction strategies for the differ- 
ent types of leptoquarks. The 'stransverse' mass variable which appears in the table, 
Mt2, has been defined previously in [17], for the case of identical semi- invisible pair 

decays as: 

M T2 = min {max (M T ,M^)} , (3.1) 

where the minimisation is taken over ^ T and fi T , the transverse momenta of the 
invisible particles, with the constraint that their sum equals p/r, the total missing 
transverse momentum, and Mt and M' T are the transverse masses calculated for the 
two decay chains. We assume that the invisible particles are massless and use the jet 
masses in our definitions of Mt2- The new variables M^ n and M min will be defined 
in section 13.4.1. 



We present our analysis of the mass reconstruction techniques for each pair- 
production decay mode separately, initially at parton level and then at detector 
level, including discussion of the relevant backgrounds. We focus on the S singlet, 
Si triplet and S1/2 doublet and outline how to generalise the strategy to all the 
leptoquark multiplets. 

It is evident from Tables [| and |] that the leptoquark decay widths are generally 
much smaller than the resolution of the detector components, and hence our analysis 
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is not sensitive to the decay widths. Throughout what follows we have set the 
leptoquark couplings to fermions to the value g = 0.1. This value represents an 
estimate of the leptoquark couplings to third generation quarks and leptons, derived 
using the measured fermion masses as explained in [6]. The resulting width-to-mass 
ratio for the leptoquarks corresponding to this coupling, according to Eq. ( |2.10| ), is 
£>(10- 4 ). 

We use the HERWIG++ implementation of the model to generate a number of 
events corresponding to an integrated luminosity of 10 fb _1 of the relevant signal and 
it background samples. Subsequently we use the Delphes framework [18] to simulate 
the detector effects and assess the feasibility of reconstruction in an experimental sit- 
uation. 6 Delphes includes the most crucial experimental features: geometry of the 
central detector, the effect of the magnetic field on the tracks, reconstruction of pho- 
tons, leptons, 6-jets, r-jets and missing transverse energy. It contains simplifications 
such as idealised geometry, no cracks and no dead material. We use the default 
parameter settings in the Delphes package that correspond to the ATLAS detector. 
Critical features of our analysis are both b- and r-tagging of jets and we caution the 
reader to take into consideration that the relevant efficiencies contain a degree of un- 
certainty at this stage of the LHC experiment and for the near future. The 6-tagging 
present in the Delphes framework assumes an efficiency of 40% if the jet has a parent 
6-quark, 10% if the jet has a parent c-quark and 1% if the jet is light (i.e. originating 
from u, d, s or g). The identification of hadronic r-jets is consistent with the one 
applied in a full detector simulation. The resulting efficiencies for hadronic r-jets 
are in satisfactory agreement with those assumed by ATLAS and CMS. See [18] for 
further details. 

Throughout the analysis we apply transverse momentum cuts of at least 30 GeV. 
Since we are always working with high-transverse momentum objects, we can assume 
that pile-up arising due to secondary proton-proton collisions is under experimental 
control. See, for example, the ATLAS itH(—t bb) study in [19]. 

3.2 (tT)(tr) decay mode 

We examine the possibility of full reconstruction of the topology shown in Figure Jl], 
where we have, for example. S (S ) — > bjjjiVi and S (Sq) —> blv^fe^-, where z/ x and 
z/ 2 represent one or more neutrinos coming from the r decays and t can be either 
a muon or an electron. We can assume that the neutrinos associated with the 
decays of the rs are collinear with the direction of the jets jx,i associated with them. 
The validity of this assumption has been confirmed using HERWIG++, for leptoquarks 
of mass 1, 0.4 and 0.25 TeV, as may be seen in Figure ^|, which shows the distribution 
of SR = ■sj 5t] 2 + 5(j) 2 between the momenta of the r jet partons and the r invisibles. 
The assumption is employed in our reconstruction of any leptoquark decay mode 
containing a r-jet. 

6 Delphes is a framework for fast simulation of a general-purpose collider experiment. 
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modes 


types 


technique 


MM 


c c(0) 
Do, 1 


j T | u T , mass constraints 
=>• edge reconstruction (M^; 1 , M min , M T2 ) 


(tr)(tr) 


Do, 1 , 

q(+) CI 

°l/2i °0 


two j r | v T , mass constraints 
=>- full reconstruction 


MM 


q qyJ) 
Do, 0± , 

q(-) q'(-) 
°l/2 ' °1 


M T2 


MM 


c(+) c(") 
°1 ' °l/2 

c q'\+> 

DO, D 1 ^ 2 , 

g/(0) 


two j r || z-Vj mass constraints 
=4> full reconstruction 


MM 


°1 ) °l/2 

c/ c'(0) 
d , d x , 

c'(-) 
D l/2 


M T2 


MM 


q(-) q' 

°l/2i J 
g'(0) 


j T | u T , mass constraints 
=>■ edge reconstruction (M^, M min , M T2 ) 



Table 7: The table outlines the general reconstruction strategy for leptoquark pair- 
production for the different types of leptoquarks. For variable definitions and further 
details see the respective sections. 
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Figure 1: Pair-production of Sq leptoquarks with decay to (tr)(tr), followed by one 
hadronic and one semi-leptonic top decay. 
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Figure 2: The distribution of the distance in i?-space (5R = s/ 8rf + 6(f) 2 ) between the 
momenta of the r jet and the r invisibles in Sq pair-production for Mlq = 1, 0.4, 0.25 TeV. 



The top quark branching ratios are ~ 0.216 for the semi-leptonic e,/i modes 
and ~ 0.676 for the hadronic top modes, this is multiplied by two (since we have 
either the t or t in each of these) resulting in a ~ 0.292 factor just for the top decay 
modes. The branching ratios and cross sections for So So production depend on the 
leptoquark mass and coupling and are shown in Table |8|, where the last column is 
the resulting cross section for the whole process under study. We focus on 400 GeV 
leptoquarks since these are clearly not excluded by direct searches and still provide 
a sufficient number of events to be potentially discovered at a reasonable luminosity 
(10 fb" 1 ) at 14 TeV. 

3.2.1 Kinematic reconstruction 

The final states of S S —> ir + tr~ processes contain many decay products including 
neutrinos. If the system has a large enough number of kinematical constraints, 
such as mass-shell conditions and balance of the total transverse momentum, we 
can completely reconstruct the kinematics of the system. The numbers of unknown 
variables and constraints are summarised in Table ^ for each decay pattern of the 
tops: (1) both tops decay hadronically, (2) one top decays semi-leptonically and 
another hadronically and (3) both tops decay semi-leptonically. As mentioned above, 
we assume r neutrinos are collinear to the r-jets, leaving two unknown parameters 
associated with the taus, namely the energy ratios Z{ (i = 1,2, > 1) which are 
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M So (GeV) 


a(pp -)> Sb-So) (pb) 


BR(tr) 


aitrtr — > bbjj£(= e, h)vtt) (pb) 


174.2 (= m top ) 


141(1) 


0. 


0. 


250 


24.3(3) 


0.34 


0.729 


400 


2.000(7) 


0.567 


0.188 


500 


0.561(6) 


0.606 


0.06 


1000 


5.94(7) xlO- 3 


0.65 


7.3 x 10" 4 


M top (GeV) 


a(pp — > tt) (pb) 




a(ti 66jj^(= e,/i)z/) (pb) 


174.2 


834(1) 




242 



Table 8: SqSq total cross section at the LHC at 14 TeV pp c.o.m energy, branching 
ratio to tr and remaining cross section taking into account the top branching ratios. The 
corresponding tt values are given for comparison. 



Decay type 


# of unknowns 


# of constraints 


(1) had, had 


1 + (0 + 2)iV 


(2 + 2)V 


(2) had,lep 


1 + (4 + 2)iV 


(5 + 2)AT 


(3) lep,lep 


1 + (8 + 2)iV 


(8 + 2)iV 



Table 9: The numbers of unknown variables ([rang], \y from top], [energy fraction of tau]) 
and constraints ([mass-shell conditions], [balance of missing momentum]) in N events of 
each decay type. The mass-shell conditions that could constrain the unknown variables 
are counted here, i.e. the mass-shell conditions on So, leptonic top, W and v from leptonic 
top decay. 

denned (neglecting masses) by: 

Pn = ZiPji 

Vu x = P n ~ Pn = & ~ l)p k , (3.2) 

where p n , Pj. and p Vi are the four- momenta of the r leptons, r-jets and r neutrinos, 
respectively. With this assumption, the unknown variables in Table ^] are the mass 
of the leptoquark, the four momenta of neutrinos from leptonic top decays and the 
energy fractions associate with the neutrinos from the tau decays. The mass-shell 
conditions that could constrain the unknown variables are counted in Table || i.e. the 
mass-shell conditions on So, leptonic top, W and v from leptonic top decay. 

As can be seen, we can wholly reconstruct the kinematics of a single event only 
in decay types (1) and (2). In decays of type (1), it would be difficult to reconstruct 
both hadronic tops because of the large combinatorial background. Thus, we focus 
on decay type (2) and attempt to determine the leptoquark mass. As we show in 
Appendix 0, in this case one obtains a quartic equation for the energy ratio Z2, 
and hence in general up to four solutions for the leptoquark mass, at least one of 
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which should be close to the true value if the visible momenta and missing transverse 
momenta are well measured. 

3.2.2 Parton-level reconstruction 

We first perform the (tr)(tr) analysis of the hard process (no initial or final state 
radiation, no underlying event) at parton-level without considering experimental or 
combinatoric effects, to examine its feasibility. For the majority of cases there are 
only two physical, approximately degenerate, solutions, which are close to the true 
leptoquark mass. The numerical solution of the quartic equation sometimes fails to 
yield roots. The results for true leptoquark masses M So = (0.25, 0.4, 1.0) TeV are 
shown in Figure |[ which includes histograms of the solutions obtained for 10 3 events. 
The histogram includes a bin at where the events without solution are placed. These 
amount to about 10% of the total events. At this level the reconstruction technique 
provides a good estimate of the leptoquark mass for all the trial true masses, lying 
within a few GeV of the true mass. 




500 1000 1500 2000 500 1000 1500 2000 500 1000 1500 2000 



Mm (GeV) Mm (GeV) Mlq (GeV) 

Figure 3: Histograms of the solutions obtained at parton level for So^o and decay to a 
semi-leptonic top, a hadronic top and two tau leptons, for Ms = (0.25, 0.4, 1.0) TeV (from 
left to right respectively). The zeroth bin contains the events where no real solution has 
been found. 

3.2.3 Experimental reconstruction 

We consider an Sq leptoquark with mass Ms = 400 GeV, for which the cross section 
for production and decay into the topology of Figure [l|, S S — > trtf — > bbjjiuTf, 
is a = 0.187 pb. The most significant background in this scenario is tt production, 
with two extra jets misidentified as rs and subsequent decay of the tops into bbjjiu. 
The cross section for this process is 242.4 pb, overwhelming to begin with. There is 
also potentially an irreducible Htt — > rfti background which, for a Higgs of mass 
Mn = 115 GeV, has a cross section of approximately 65 fb. Since one of the main 
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rejection mechanisms is the reconstruction through the solution of the kinematic 
equations, we do not expect this background to contribute significantly. 

We simulate the events with QCD initial state radiation (ISR), final state radi- 
ation (FSR) and underlying event (UE). We use the default jet algorithm provided 
by the Delphes package for the ATLAS configuration, anti-kt with R = 0.7. We 
then demand a set of relatively loose cuts on the full it and SqSq samples, since in 
a real experiment we wouldn't be able to separate the different decay modes of the 
top quark or So leptoquark. The cuts applied are the following: 

• the existence of a lepton in the event, being either a muon or electron, with 
PT >e > 30 GeV. 

• a minimum of 6 jets. 

• the missing transverse momentum in the event, B}t > 20 GeV. 

• two r-tagged jets, with the extra requirement that they both have pt, t > 
30 GeV. 

• no jets tagged as both b- and r-jets simultaneously. 

We also require that the highest-p^ lepton is at a distance SR > 0.1 from the r- 
tagged jets, since electrons may create a candidate in the jet collection as well as the 
lepton collection. The analysis then breaks up into different branches according to 
the number of 6-tagged jets in an event. 

• two b-tagged jets: we look for one or two further jets (with > 30 GeV) that 
form an invariant mass close to the top mass, within 20 GeV. One 6-jet is then 
associated with the semi-leptonic top decay and the other with the hadronic 
top decay. 

• one or no b-tagged jets: when there is one 6-tagged jet we check whether it 
will satisfy the top mass conditions with any other (one or two) remaining jets, 
otherwise we associate it with the semi-leptonic top. If so, we look for any two 
or three jets that satisfy the top mass conditions, and form the hadronic top 
within a 20 GeV mass window. For the remaining fe-jets (or if there are no 
6-jets) we look for the remaining highest-pr jets. Any jets that are found in 
this way and called 6-jets are required to have a px,b > 30 GeV. 

No solutions are found in the sample of 70 signal events passing the cuts, if we require 
the ratios Zj to be purely real. Hence, the solutions to the quartic equation for the 
momentum ratios Zi, described in Section 3.2. lj , are now allowed to be complex in 



order to provide some signal, since even true leptoquark events are smeared and 
distorted by detector and QCD effects. We use the real part of z 2 as an input to 
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the calculation of the rest of the kinematic variables. This is reasonable since the 
experimental effects are expected to 'smear' the position of the true value of z 2 in 
all directions in the complex plane. The effect is shown in Fig. |], where we plot the 
real and imaginary parts of z 2 for the events that have passed the kinematic cuts. 
Evidently, there is a concentration of true solutions around the positive real axis, 
an effect exemplified by Fig. |5], where we show the ratio of the real part of z 2 and 
its modulus. We have further demanded that the resulting momentum fractions are 
physical: TZ(zi i2 ) > 1, resulting in only real solutions for Ms Q . Figure || shows a 
reconstruction plot for leptoquarks of mass M So = 400 GeV. Note that each event 
was given weight 1, distributed evenly amongst the solutions it yields. In the case of 
complex z 2 , we assume there is one solution corresponding to the complex conjugate 
pair. 
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Figure 4: The plot shows the complex values of the solutions for z 2 after solving the 
quartic equation for the events that have passed the experimental cuts. There exists a 
higher concentration of events about the positive real axis. The number of entries is 280 
(4 solutions included for each of the 70 events). 

Although the cuts applied are relatively weak, most of the rejection comes from 
the requirement of two r-tagged jets. The background does not produce solutions in 
the physical region often enough to be significant. 

3.3 (qis)(qv) decay modes 

We can obtain the mass of the leptoquarks when both of them decay into bv or tu 
using the M T2 variable (Eq. (|3.1|)). Examples of these decay mode are S0S0 — » bubu 
and 5 , {/2^'i/2 ~~ ^ t^tv. 

3.3.1 Parton-level reconstruction 

At parton level, the (tv)(tv) and {bv){bv) decay modes are similar and hence we 
consider only the latter here. We first construct the M^ 2 variable using the parton- 
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Figure 5: The plot shows the ratio of the real part of zi and its modulus. The peak 
close to 1 demonstrates the clustering of the real positive solutions about the real axis and 
justifies the use of the real part as an input to the rest of the calculation. The number of 
entries is 280 (4 solutions included for each of the 70 events). 
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Figure 6: Experimental reconstruction of the SqSq — > trif — > bbjj£(= e, \i)vtt mode 
using the method described in the text. Note that each event has weight 1, distributed 
evenly amongst the solutions it yields. The signal is shown in red (35 entries) and the ti 
background in blue dashes (3 entries). 
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level 6-quark 4-momenta, in the absence of any experimental effects, ISR or FSR. 
The result is shown in Figure [7] for Mlq = (0.25, 0.4, 1) TeV, confirming the expected 
sharp edge in these idealised conditions. 
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Figure 7: The parton-level Mti distribution constructed for the SqSq —> bvbv using the 
true 6-quark momenta, for Mlq = 1,0.4,0.25 TeV. 



3.3.2 Experimental reconstruction 

As before, we use the Delphes framework to simulate the detector effects, with the 
settings stated in Section [TTJ We demand two 6-tagged jets in both the q = b and 
q = t cases. In the latter we search for combinations of 1 or 2 jets with the 6-tagged 
jets which form the top mass within a window of 10 GeV. We now require the 
following cuts for the (bv)(bv) case, on the full S S sample: 

• two b-tagged jets with pr,b > 120 GeV each. 

• no electrons or muons in the event. 

• Missing transverse energy E}t > 250 GeV. 

For the {tv){tv) case we require the following cuts on the Sy^Sy^ sample: 

• two b-tagged jets with p Tib > 80 GeV each. 

• no electrons or muons in the event. 
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• Missing transverse energy B}t > 260 GeV. 

The resulting My 2 distributions for the signal (blue) and ti background (red) can be 
seen in Figure || The (tv) {tv) mode appears to be more challenging to reconstruct 




Figure 8: Experimental reconstruction of the SqSq — > bv T bi7 T mode (left, 121 background 
events, 125 signal events) and Sy 2 Sy 2 — > tvtv mode (right, 39 background events, 48 
signal events) using Mt2- ISR and FSR and the underlying event have been included in 
the simulation. The signal is given in red and the ti background in blue dashes. 

than the (bu)(bu) mode. This is due to the fact that the ti background is very similar 
to the signal and the difficulties that are presented in reconstructing hadronic tops. 
Nevertheless, as the results show, it may be possible to observe an excess over the 
Mt2 distribution of the background and provide an estimate of the mass. 

3.4 (q'T)(qv) decay modes 

One possible event topology for the S So — > bvtr~ processes is shown in Figure |9[ 
When the top decays hadronically, the system has two neutrinos, one from an So 
decay and another from a r decay. If we can reconstruct the hadronic top correctly, 
we can simply use Mt2 to obtain the mass of the leptoquarks. Similar topology is 
present in the S[^S[^ ~~ (^)(^) decay mode. 

It is known that the information from Mt2 is the same as that from the 'minimal 
kinetic constraints', in events where two identical particles decay to missing particles 
with the same mass [20-22]. As discussed in Section in this type of event, we 
can take advantage of the fact that, to a good approximation, the neutrino from a 
r decay is travelling almost collinearly to the r jet in the lab frame. By including 
this constraint, we can define kinematical variables, M min and M^, which perform 
better than M^ 2 at parton-level, as we will show in the following subsections. 
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Figure 9: Pair production of So leptoquarks with decay to (tT)(bu), followed by hadronic 
top decay. 



3.4.1 Kinematic reconstruction 

In the above approximation we can write (neglecting masses) 

Vv-r = w Pj (w > 0) . 



(3.3) 



The second neutrino comes directly from the So decay associated with a 6-jet. The 
transverse components of the momentum of this neutrino are constrained by 



Pmiss - Wp 



3 ■ 



(3.4) 



There are two unknown parameters left, w and p z v . In terms of these, we define two 
invariant mass variables: 



and 



m 2 T (w) = (p t + (1 + w)pj) 2 = m 2 t + 2(1 + w)p t ■ Pj 



m 2 bu (w,p z u ) = (p b + p. 



(3.5) 



= 1E h \j (p miss - wpj) 2 + (p z v ) 2 - 2p b ■ (p miss - wpj) - 2p z h pl . (3.6) 

Note that m tT does not depend on pf, and is a monotonically increasing function of 
w because p t ■ Pj > 0. We can now define two M<r 2 Tike variables: 



M min = min[max{m 4r ,m 6i ,}] > M T2 



and 



M, 



X= min [m hu ] , 

m t T=m bv 



(3.7) 
(3.8) 
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Figure 10: Parton level distributions of (left), M m - m (centre) and Mti (right) for 

(bu)(tr) (solid curve) and (tv){br) (dashed curve). 



where minimisation is taken for all possible (w, pi). By construction, both these 
quantities have an upper bound equal to the leptoquark mass: 

M So > M min , M So > M^l . (3.9) 

Furthermore we show in Appendix [E| that 

M™>M min . (3.10) 

3.4.2 Parton-level reconstruction 



Figure |i~0l shows the parton-level distributions of M ^ , M min and M T2 . Here we 
generated 1000 events and took only the true combination of the jet assignment. As 
can be seen, all the distributions have clear edge structures at the input leptoquark 
mass of 400 GeV. 

In order to compare these variables we took the differences, shown in Figure pTT]. 



The relation Af£& > M min > M T2 is seen to hold on an event-by-event basis. This 

fbal 
'min 



implies that M^ n and M min are more powerful than Mj>2 for determining the mass 
of the leptoquark, at least at parton level. 



3.4.3 Experimental reconstruction 

The settings for experimental reconstruction used for the Delphes fast simulation 
remain unaltered in the present analysis (see Section |3.1| ). We apply the following 
event selection cuts to the full S S signal and the tt background: 

• at least four jets found in each event. 

• exactly one r-tagged jet with pr > 120 GeV. 

• no, one or two 6-tagged jets with p? > 60 GeV. 

• missing transverse energy, BJt > 200 GeV. 
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Figure 11: Parton level distributions of - M min (left) and M min - M T2 (right) for 

{bv){tr) (solid curve) and (tz/)(6r) (dashed curve) 



For the fe-jet originating from the leptoquark decay, we choose the highest-pr fe-tagged 
jet when there are two 6-tagged jets and the highest-pr jet (excluding the r-tagged 
jet) when there are no ^-tagged jets. We use all the remaining jets with p? > 30 GeV, 
(not identified as the 6-jet from the leptoquark) to search for one, two or three jets 
that form an invariant mass close to the top mass, within a 20 GeV window. We 
apply the additional constraint that the difference between the pt of the r-tagged 
jet and the pt of the 6-tagged jet, px, T — PT,b > — 10 GeV. This eliminates a high 
fraction of the ti background since the rs in that sample originate from the W decay 
and are expected to have lower p? on average than the bs that originate directly from 
the top. On the contrary, in the leptoquark signal the r and b transverse momenta 
are expected to be more equal on average. 



The resulting distributions are shown in Figure |12| Due to the low number of 
events passing the selection cuts, it is not obvious whether the observable per- 
forms better than M m j n and Mt2- We checked, however, that the three distributions 



satisfy the same inequalities presented in Figure |TT] for the parton-level reconstruc- 
tion. 

The reconstruction strategy for the (tu) (br) mode follows the technique described 



in Section |3.4.1| for the (bu)(tr) case, with the simple replacement b -h- t. The 



assignments of 6-jets and top-jets is performed in the same way as in the [bv){tr) 
analysis, with the following cuts applied to the full S\_2s\_2 sample: 

• at least four jets found in each event. 

• exactly one r-tagged jet with pr > 190 GeV. 

• no, one or two 6-tagged jets with px > 40 GeV. 

• missing transverse energy, E}? > 120 GeV. 
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Figure 12: Reconstructed distributions M^ n (left), M m [ n (centre) and Mt2 (right) for 
(bu)(tr) signal (red) and ti background events (blue dashes) including ISR, FSR and the 
underlying event. There are 37 signal events and 19 background events in all plots. 



There is also a cut on the reconstructed hadronic top jet, of pr > 120 GeV and 
that its invariant mass lies within 20 GeV of the top mass. The results are shown 
in Figure [13]. Note that the background that would be present due to the Sy£ 
leptoquark has not been included. 

Although at parton- level, the variable performs better than M min and M^ 2 , 
it seems to become unstable after including experimental errors, with some events 
failing to produce a value within the range of the plots shown in Figure O. The 
origin of the instability is the additional assumption of the leptoquark masses being 
equal, which is satisfied at parton-level (up to small width effects) but does not hold 
exactly after detector simulation. For the events for which no solution is found, we 
assign = M m i n . Even after this readjustment, there are a few events for which 
a solution for M^, is found and lies outside the region shown. Therefore, Mt2 and 
A^min appear to be preferable as experimental observables. 

3.5 (6r)(br) decay mode 

3.5.1 Kinematic reconstruction 

The (6r)(fer) mode can be fully reconstructed if one again assumes collinearity of 
the r-jets and r-neutrinos: p T ^ = Zip^i (i = 1,2, Zi > 1). This implies that the 
missing momentum from each r can be written as pji = (zi — l)pj,i- Hence, we may 
write the following equalities for the components of the measured missing transverse 
momentum: 

Pmiss = p X n (Zl ~ 1) + P X j2 (z 2 - 1) 

Pmiss = Pjl( Z l ~ 1) + Pj2^2 - 1) • (3.11) 
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Figure 13: Reconstructed distributions M^?* (left), M m i n (centre) and Mt2 (right) for 
the (tv){br) signal (red) and ti background events (blue dashes) including ISR, FSR and 
the underlying event. There are 68, 72 and 72 signal events (left to right) and 22, 23, 23 
background events (left to right). 



The above equations may be written in matrix form and inverted to give: 

Vj 2 "miss "_/'2.rmiss 

Z l J- H Z y V Z ) 

PjlPj2 - I'fll'fl 



7 1 "miss V jl" miss / r>\ 

z 2 = 1 - y y—^— ■ 3 - 12 

PjlPj2 - I'/ll'fl 



Now the invariant mass of each of the two leptoquarks may be written as m| 
(pb +Pt) 2 , resulting in the following expression: 



s 



2z i p bi -p ji , (3.13) 



where we have neglected the r and 6-quark mass terms. Using Eqs. ( |3.12 ), we obtain 



two values of ms on an event-by-event basis. At parton level, with the correct jet 
assignments, these solutions approximate the leptoquark mass very closely, up to the 
collinearity approximation. 



3.5.2 Experimental reconstruction 

As before, the Delphes framework has been used with identical settings. The fol- 
lowing cuts have been applied to the -> (6f)(6r) mode: 

• at least 4 jets present in the event. 

• two r-tagged jets with pt > 140 GeV. 

• missing transverse energy B}t > 140 GeV. 
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We accept events with no, one or two 6-tagged jets. If there are less than two 6-jets, 
we search for the highest-pr non-tagged jet(s) to obtain two 6-jets. We apply a cut 
of pt > 50 GeV on these. There are two possible assignments of the br combination, 
resulting in a total of four solutions. The resulting distribution for the mass solutions, 
as described in Section |3.5.1|, is show in Figure O. The tt background appears to be 




Figure 14: Experimental reconstruction of the (6r)(6r) mode using the method described 
in the main text. ISR and FSR have been included in the simulation. Note that each 
solution has weight 0.25. The signal is shown in red (18.75 entries) and the ti background 
in blue dashes (1.75 entries). 



under control, with 1.75 entries in the mass histogram. 

We also considered the bbjj background, for which we generated events using 
Alpgen v2.13 [23], applying the px cuts on the four partondevel objects. We con- 
cluded that we can safely ignore this background since the expected number of events 
with two r-tagged jets was O(l), before applying any restrictions on the missing 
transverse energy, E}t- Note that the backgrounds to this decay channel originating 
from the other members of the leptoquark multiplet have not been included. 

3.6 Determination of quantum numbers 

In the ideal scenario where all of the decay modes of a leptoquark multiplet are seen, 
the quantum numbers can be deduced without ambiguity. For example if we only 
observe combinations of (tr) and (bu) decay modes, then the only likely candidate 
is an 5*o singlet. If in conjunction with these decay modes we observe (br) and 
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(tv) decay modes, with corresponding total rates, then we might guess that we have 
observed the S\ multiplet. 

The issue is more complicated if some decay modes are missed. For example if 
only the (tr) (tr) decay mode has been seen, we might assume that we have observed 
the pair production of an S' leptoquark. However, we might have observed the 
(tr)(tr) decay of an S^t? leptoquark pair and missed the more challenging (tv)(tv) 
mode of the S'^ leptoquarks. In this case we would need to examine the helicities 

and charges of the decay products: the decays to i L f L and i R f R whereas the S' 
decays to t^fi and tijR. Since we can reconstruct all decay products of the top and r 
without combinatorial ambiguity, using measured leptoquark mass as an input, there 
is hope that we could measure top [24-26] and r [27,28] polarisations simultaneously. 
This would allow us to distinguish these two cases. We leave investigation of the 
feasibility of this to future work. 

4. Conclusions 

If strongly-coupled dynamics solves the hierarchy problem of electroweak symmetry 
breaking, the question arises of how best to discover it at the LHC Existing con- 
straints coming from electroweak precision tests tell us that, at least at low energies, 
any model of strong dynamics must be a lot like the Standard Model (with perhaps 
a faint hope of observable deviations in the Higgs sector [29,30]). 

Existing constraints coming from flavour physics are somewhat different in that, 
while the data are certainly consistent with the Standard Model, naturalness argu- 
ments suggest that strongly-coupled theories should differ from the Standard Model 
in the flavour sector. Indeed, fermion masses should arise via mixing between ele- 
mentary and composite fermions of the strongly-coupled sector. 

If that is so, then composite leptoquarks (or diquarks) may also appear, cou- 
pled predominantly to third-generation fermions. These would provide a spectacular 
signature at the LHC. Their Standard Model quantum numbers imply that they 
would be produced strongly as conjugate particle-antiparticle pairs, decaying into 
third-generation quarks and leptons in the combinations summarised in Table [7[ 
We have proposed a number of new experimental search strategies adapted to these 
characteristic final states, also summarised in Table [7], and implemented the relevant 
processes in the HERWIG++ event generator in order to study their effectiveness in 
the presence of QCD radiation, backgrounds and the underlying event. We used the 
Delphes detector simulation to assess the effects of b and r tagging efficiencies and 
detector resolution. For definiteness we assumed a leptoquark mass of 400 GeV and 
an integrated pp luminosity of 10 fb _1 at 14 TeV. 

In the case of decays of leptoquark pairs to (qT)(qr) where q = t or b, the 
approximate collinearity of the missing neutrinos and jets from the tau decays allows 
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full reconstruction of the leptoquark mass, even when one top decay is semileptonic. 
In the former case there is a quartic ambiguity in the resulting mass, although not all 
of the solutions are real. After detector resolution smearing, the correct solutions for 
the momentum fraction z 2 may be complex, but we found that using the real parts 
provides a fair estimate of the mass, with resolution of the order of ±150 GeV. For 
(6r)(6r) the only ambiguity is combinatoric but the mass resolution is similar. In 
both cases the expected background from qqjj is small after cuts and reconstruction. 

For decays to (tr)(bv) or (tu)(br), we have proposed an edge reconstruction strat- 
egy similar to those developed for supersymmetric models, but using mass variables 
M^l and M min that are in principle superior to the classic 'stransverse mass' My 2 . 
However, given the limited statistics expected, the difference in performance between 
these variables was not obvious. We found cuts to reduce the background from tt to 
manageable levels, but the edge reconstruction remains challenging without higher 
statistics. For (qv)(qv) the story is similar for edge reconstruction in M T2 , the case of 
q — t being the more difficult owing to the similarity of the distributions of the signal 
and ti background. But even in that case a clear excess over background should be 
visible and would give a rough estimate of the leptoquark mass. 

In the event that a discovery is made, one might ask to what extent this provides 
proof that electroweak symmetry breaking is driven by strongly-coupled, composite 
dynamics. After all, one can easily imagine weakly-coupled theories with such states, 
for example, third generation squarks in R-parity-violating supersymmetric models. 
Ultimately, TeV-scale compositeness can only be revealed by experiments probing 
significantly higher scales; for that, we shall have to wait some time. In the meantime, 
the discovery of leptoquarks coupled to third generation fermions and their de facto 
consistency with the multitude of existing flavour experiments would imply very 
strong bounds on the couplings to other fermions. The scenario in which the observed 
fermions are partially elementary and partially composite provides, as far as we 
know, the only mechanism in which the required suppression can be automatically 
achieved. Moreover, it gives a prediction for the size of the other couplings, some 
of which are not far from current bounds, which may then be targeted in ongoing 
flavour experiments. Though circumstantial, this would seem to be the best possible 
evidence for compositeness that one might hope for in the LHC era. 
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A. Feynman rules and diagrams 



The Feynman rules [15] relevant to the leptoquark pair-production diagrams are 
given in Figure |T5| and Figure IB. The relevant parton- level Feynman diagrams are 
shown in Figure [17] and Figure pl| for gluon-gluon and quark- ant iquark initial states 
respectively 



ySS,aij _ g s (t a yj(k 2 ~ h) t 



Figure 15: Feynman rule for the vertex scalar leptoquark-scalar anti-leptoquark-gluon. 
All momenta are incoming and arrows indicate colour flow. 
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Figure 16: Feynman rule for the vertex scalar leptoquark-scalar anti-leptoquark- gluon- 
gluon. All momenta are incoming. 

Since the couplings to light generations are suppressed in the kind of models 
we are considering, leptoquark single-production in hadron colliders can proceed 
only via 6-quark gluon fusion, as shown in Figure |19[ However this is also heavily 
suppressed due to the low 6-quark PDF and the small couplings to fermions, and can 
be neglected. 



B. Cross sections 



Using the given Feynman rules, the cross section for scalar leptoquark pair-production 
is completely determined, with the only unknown parameter being their mass. The 
partonic differential cross sections with respect to the leptoquark scattering angle in 
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Figure 17: Feynman diagrams relevant to scalar leptoquark pair-production with gluon- 
gluon initial states. 
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Figure 18: Feynman diagram relevant to scalar leptoquark pair-production with guark- 
antiquark initial states. 
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Figure 19: Feynman diagrams relevant to scalar leptoquark single production. Solid lines 
with an arrow indicate quarks, lines without an arrow indicate leptons. 



the partonic centre of mass frame, 9 are given by: 
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where the 'boost factor', 0, is given by (3 = yj 1 — 4M|q/s, with s the square of the 
parton-parton centre of mass energy. The integrated partonic cross sections are then 
given by: 

1 + 



na 2 



*S = ggf < /?(41 - 31/3 2 ) - (17 - 18/3 2 + /3 4 ) log 



P 



?irry 2 

*h = • (R2) 

Note that the production cross sections at leading order in QCD are exactly equiva- 
lent to the supersymmetric top partner pair-production cross sections. This fact was 
verified directly using the implemented supersymmetric model in HERWIG++. 

We obtain the hadronic cross sections in the usual way, i.e. by folding with the 
parton density functions and integrating with respect to the momentum fractions 

a ss = f dx 1 dx 2 f hl ,i( x i'S)fh2,j( x 2iS)°ss& ' ( B - 3 ) 

where Eh lt h 2 is the hadronic centre-of-mass energy, s = x\XiE\ h , and fh a ,i( x k, s) is 
the distribution function for parton i in hadron h a , evaluated at scale s and momen- 
tum fraction x^. 

The vertices appearing in Figure ^ and Figure |15] have been implemented in 
the HERWIG++ event generator [10, 11] and the diagrams contributing to the scalar 
leptoquark pair-production are reproduced automatically therein. This yields a cal- 
culation for the cross section which, as a check on the implementation, we compare 



in Figure ^0] with direct integration of Eq. ( |B.3| ) using the adaptive integration al- 
gorithm VEGAS [31]. For this comparison we used the MRSTMcal (also known as 
LO**) parton distribution functions [32] in both cases. The slight discrepancy is due 
to the difference between the internal definition of the strong coupling constant in 
HERWIG++ and the one associated with the MRSTMcal PDFs. 



C. The effective Lagrangian for decays of derivatively- coup led 
lepto quarks 

The Lagrangian for derivatively-coupled conjugate fields, which appears in Eq. (|2.5| ), 
also contains terms involving the conjugate fields, such as: 

£§> 1/2 - i C Rl»T LP ^S'$ . (C.l) 

To manipulate the above expression for the case of on-shell S[, 2 decays as we did in 
Eq. (|2.7| ), we need to show that: 

*W = m^l R , (C.2) 
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Figure 20: The variation of cross section for scalar leptoquark pair-production at the 
LHC (14 TeV pp centre-of-mass energy) using HERWIG++ and direct integration using the 
VEGAS algorithm (MRSTMCal parton density functions). 



where \I/ is a 4-component spinor and R = (Pl^^/) c . This can be done by using 
the following identities [33]: 



c-\ = -ilc- 1 



(C.3) 



and hence: 



So the necessary effective Lagrangian for the decay is given by: 



(C.4) 



C eff ~ m t t c L T L S'$ 



(C.5) 



as required. The full list of effective Lagrangians for the primed leptoquarks, from 



- 32 - 



which the decay modes and couplings in Tables |2| and [| can be derived, is given by: 
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(CIO) 



where we have defined: = (Si =F iS'i )/V2 (and equivalent definitions for 

S'^) and S'^ = S'i 3 \ We have also used the fact that the doublet leptoquarks 



may be written as a vector S[/ 2 = (S'^ 2 \ S'^ 2 ). We have set the quark and lepton 
couplings to equal, g q = g e , 7 however these can be reinstated trivially by replacing 
g — >■ g q where a quark mass term appears and g — )■ g where a lepton mass term 
appears. 

Note that terms appearing in this Lagrangian are no longer SU(2)l x Uy gauge- 
invariant. This is consistent since these terms would appear in the Lagrangian after 
electroweak symmetry breaking and vanish as the fermion masses tend to zero. The 
Lagrangian is, of course, U(l) em gauge- invariant. 



D. (tr){tr) reconstruction method 

In terms of the momentum ratios Zj defined in Eq. (|3.2|), the conditions for balancing 



7 The implementation in the upcoming HERWIG++ version 2.5.0 also includes this simplification. 
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the total missing transverse momentum can be written as 



v\v x n -v x jx vl x 
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The mass-shell conditions, except for p^ =0, can be written as 
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(D.l) 
(D.2) 

(D.3) 
(D.4) 
(D.5) 
(D.6) 



where m t is the reconstructed mass of the hadronic top and m t is the assumed mass 
of the semi-leptonic top. By eliminating z\ and ms from eqs. ( P-l|) , ( P-5|) and (p.6| ), 
one obtains 



, Pt ' Pjl X . | 
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Using a vector p 
as 

where 



u 3 = -(Pb+Pi) - Ph 
{E n p* p\ plX eqs. Q, O, O and ((T^D can be recast 



AP 



/ E 



A 



1 



-Pt 
-Pb 



-Pi 

V 

-Pb 



-Pt \ 
-Pb 



z 2 E h -z 2 p% + (j>fPh)/P X h ~ z ^Ph ~ z ^P Z h 



\ 



V 

Ph 



-Ph 







/ 



and 



2 ' 



£4 and m 4 are defined as 

U 
U4 



mt —mt—mt 



-Pb-Pu t 3 + u 3 z 2 , t 4 + U 4 Z 2 



(Pmiss +P%)Ph ~ (Pmiss + P? 3 )l&, 



T V V 7' 

PhPh-PhPh- 



(D.10) 
(D.ll) 

(D.12) 

(D.13) 
(D.14) 



From Eq. ( p.lOj ), one can determine p U[ as a function of z 2 . Finally, one can determine 
z 2 from the mass-shell condition 

(A^S) 2 



0. 



(D.15) 



This provides a quartic equation for z 2 , and one can find up to four real solutions in 
the physical range z 2 > 1. One can then obtain ms by substituting z 2 into Eq. (D.6). 
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E. (q f r)(qv) reconstruction method 

Given w in Eq. (|3.6|) , one can minimise m bu (w,pl) in terms of p z v . The result is: 

= 2|p fc ||p miss - wpj\ - 2p b ■ (p miss - wpj) 

= K»F, (E.l) 



where 



Pu = i — i Pb 

IP&I 

and m^(w) is the transverse mass of the bv system. This allows us to calculate M m ; n 
by one parameter minimisation 

Mmin = vam[max{mt T (w), {w)})- (E-3) 

w 

Since m tT (w) is a monotonically increasing function of w, if m ir (0) > m^(0), 

M min = m tr (0) . (E.4) 

Furthermore, since there exists a value p z v which fulfils m tr (0) = mbj,(0,]3^), one finds 

M^ n = m tT (0) (E.5) 

If m^(0) > m tT (0), we have to search for other values of w. For the true w and p z v , 
say u>* and pi*, we have 

m£>*) < m bv {w\p z *) = m tT (w*) . (E.6) 

This assures existence of w which satisfies the relation m b ^{w) = m tT (w). By scanning 
w from to w, one finds 

MX = O) , (E.7) 

and 

M min = min [m%{w)] . (E.8) 

Hence we have: 

> M min . (E.9) 
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